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Abstract 

We prove that the matrix of capacitance in electrostatics is a 
positive-singular matrix with a non-degenerate null eigenvalue. 
We explore the physical implications of this fact, and study the 
physical meaning of the eigenvalue problem for such a matrix. 
Many properties are easily visualized by constructing a "poten- 
tial space" isomorphic to the euclidean space. The problem of 
minimizing the internal energy of a system of conductors un- 
der constraints is considered, and an equivalent capacitance for 
an arbitrary number of conductors is obtained. Moreover, some 
properties of systems of conductors in successive embedding are 
examined. Finally, we discuss some issues concerning the gauge 
invariance of the formulation. 
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1 Introduction 

The concept of capacitance and the matrix of capacitance has 
been studied from several points of view PQ-[H]. On the other 
hand, the theory of positive matrices and operators is extensively 
used in branches of Physics such as the mechanics of rigid body 
motion, quantum mechanics (2 1 1 122] . and other more advanced 
topics |16]-|20], Nevertheless, the employment of the theory of 
matrices and operators to study the matrix of capacitance is 
rather poor [231- 26 . In particular, no physical meaning is usu- 
ally given to the eigenvalue problem of the matrix of capacitance. 
We shall see that the theory of positive matrices and operators 
could provide another point of view that enlighten many math- 
ematical and physical properties of systems of electrostatic con- 
ductors. The paper is distributed as follows: section defines 
the electrostatic system of conductors that we intend to study, 
and establish the notation and properties necessary for our sub- 
sequent developments. In Sec. along with Appendix [X] the 
main goal is to prove the positivity of the matrix of capacitance. 
Sec. [4] discusses some subtleties with respect to the gauge in- 
variance of the formulation. Section [S] along with appendix iBl 
explores the physical implications of the positivity of the matrix 
of capacitance. This is done by constructing a "space of po- 
tentials" with inner product in which the matrix of capacitance 
represents an hermitian positive operator. Section [6] studies the 
problem of minimization of the internal energy of a system with 
constraints, and an equivalent capacitance is defined for a system 
with arbitrary number of conductors. On the other hand, config- 
urations of conductors that are successively embedded deserves 
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special attention because many simplifications are posible, and 
this is the topic of Sec. [7|and appendix [Cl Section[8]summarizes 
our conclusions and appendix|D]contains suggested problems for 
readers. 

2 Basic Framework 

This section summarizes some properties of the matrix of ca- 
pacitance obtained in Ref. [27]. They are the framework of 
our developments in the remaining sections. Let us consider a 
system of N conductors and an equipotential surface that sur- 
rounds them, such equipotential surface could be the cavity of 
an external conductor. The potential on each internal conductor 
is denoted by ifii, i = 1, 2, . . . , N. (see Fig. Q]). We define a set 
of surfaces Si sligthly bigger than the surfaces of the conductors 
and locally parallel to them, is an unit vector normal to the 
surface Si pointing outward with respect to the conductor. The 
potential of the equipotential surface is denoted by yjv+i and 
we define a surface Sn+i slightly smaller and locally parallel to 
the surface of the equipotential. The charges on the conductors 
are denoted by Qi with i — 1, N and if there is a cavity of 
an external conductor in the equipotential surface we denote the 
charge accumulated in such a cavity by Qjv+i, the unit vector 
njv+i points inward with respect to the equipotential surface. 
Finally, we define the total surface St = Si + . .. + Sn+i and the 
volume Vs T defined by the surface St i.e. the volume delimited 
by the external surface Sjv+i and the N internal surfaces Si. 




Figure 1 : N conductors surrounded by an equipotential sur- 
face. The volume Vs T is the region in white. 

Let us define a set of dimensionless auxiliary functions /; that 
obey Laplace's equation in the volumen Vs T with the boundary 
conditions 

V 2 /,— 0, = (i,j = l,...,N + l). (1) 

The uniqueness theorem ensures that the solution for each fj 
is unique in Vs T . The boundary conditions (fTJ indicate that 
the fj functions depend only on the geometry. Since the func- 
tions acquire constant values on the surfaces Si with i = 



2 



R. A. Diaz. W. J. Herrera 



1, . . . , N + 1, it is clear that V/j is orthogonal to these surfaces. 
The functions fj have some properties [27] 

JV+l 

J2 Si = 1; V/j (Si) • n* = (1 - 2<5 if ) ||V/j (Si)|| ; < /j < 1 

3=1 

(2) 

From these auxiliary functions we can construct a matrix that 
provides a linear relation between the set of charges Qi and the 
set of potentials tpi in the following way 



a 



£o i V/j -n z dS = e [ V/ • V/j (3) 
■/Sj ■/ v St 

JV+l 

= (4) 

3=1 

and some properties of the C?ij matrix can be derived 

JV+l JV+l 

dj = Cji, ^2 Cjj = ^2 Cij = 0, (5) 

3=1 »=1 

C« > , Cy < 0, (i^j). (6) 
The equations above are valid for i, j — 1, . . . ,N + 1. The ex- 
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and expressions for the internal electrostatic energy U of the 
system and of the reciprocity theorem can be obtained 

- JV+l JV+l JV+l JV+l 

i,j i i — 1 3=1 

where {Qi,^} and {Qi,^} are two sets of charges and poten- 
tials over the same configuration of conductors. The dj ele- 
ments constitute a real symmetric matrix of dimension (N + 1) x 
(N + 1), in which the number of degrees of freedom is N(N + 
l)/2, note that it is the same number of degrees of freedom of a 
TV x N real symmetric matrix. 

For future purposes, we shall call the matrix with elements 
dj and with i,j = 1, ... ,N the r-matrix (restricted matrix de- 
noted by C), while the CV, matrix with i, j = 1, . . . , N + 1 will 
be called the e-matrix (extended matrix denoted by C e ). 



3 Discussion of the mathematical 
properties of the matrix 

Equations (fj} and © tell us that the Cy elements are purely 
geometrical. In addition, Eqs. ([3| and ([5]) say that the e-matrix 
is a real symmetric matrix in which the sum of elements of each 
row and column is null. From Eq. © the non-diagonal elements 
of the e-matrix are non-positive. The volume integral in Eq. 
© shows that the diagonal elements Ckk are strictly positive 
for any well-behaved geometry. In particular, since Cjv+i,jv+i is 



positive, Eq. (f7c|) shows that at least one element of the form 
Ci,jv+i is different from zero (negative) for i = 1,...,jV; thus 
rewriting Eq. © in the form 



"52 di 

3 = 1 
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i,N + l 



(9) 



we see that if C^jv+i < the sum of the elements of the i— row 
of the r-matrix is positive, if C^jv+i = such a sum is null. 
Since at least one of the Ci,jv+i elements is strictly negative, 
we conclude that in the r-matrix the sum of elements on each 
row is non-negative and for at least one row the sum is positive. 
Because of the symmetry, all statements about rows are valid for 
columns. 

On the other hand, when Vs T is a connected region as in Fig. 
[1] the function fj should change progressively from its value 1 
on conductor j up to the value zero in the conductor i with- 
out taking local minima or maxima according to the proper- 
ties of Laplace's equation. According with Eq. ([2| the factor 
V/j (Si) ■ rii is positive for i ^ j and from Eq. © the non- 
diagonal dj factors must be strictly negative for a well-behaved 
geometry. This discussion is not valid when the volume Vs T is 
non-connected as in Fig. [2] we shall discuss this case in section 
[7] When dj < for i 7^ j, the discussion below Eq. (|9]), leads 
to the fact that the sum of elements in each row of the r-matrix 
is positive. 

In conclusion, for the e-matrix the sum of elements of each row 
is null. Further, if Vs T is a connected region, all matrix elements 
of the e-matrix are non-null (for a well-behaved geometry), and 
for the r-matrix the sum of elements of each row is positive. The- 
orems A and B in Appendix[X] show that under these conditions 
we find: O The e-matrix is a real singular positive matrix, its 
null eigenvalue is non-degenerate and the other eigenvalues are 
positive. © The r-matrix is a real positive-definite matrijfl. Its 
eigenvalues are all positive. © The null eigenvalue of the e-matrix 
is associated with (N + 1) —dimensional eigenvectors of the form 



4<o = if 0,^0, ■ 



, V>o) 



(10) 



4 Gauge invariance of the formula- 
tion 

We have two possible scenarios here, in the first the equipotential 
surface is the surface of the cavity of a conductor that encloses 
the others. In the second, the equipotential surface is just a 
geometrical place in the vacuum. The uniqueness theorem guar- 
antees the same solution in both cases but only in the interior of 
the equipotential surface. In the equipotential surface itself we 
can see that in the first case there is a charge Qjv+i accumulated 
in the cavity, while in the second case there is no charge in such a 
surface at all. The problem lies in the fact that the electric field 
is not well-behaved in the surface of the cavity because of the 
accumulation of surface charge [23]- [26], it is precisely because 
of this fact that we defined surfaces slightly different from the 
real surfaces on each conductor (in which V/j are well-defined). 
So all the observables (charges, potentials, electric fields) are the 
same in the interior of the equipotential surface for both scenar- 
ios, but the surface charge and the electric field differ in both 



*The non-degeneration of the null eigenvalue of the e-matrix follows 
from theorem A or alternatively from theorem B, in appendix [A] after 
establishing the positive -definite nature of the r-matrix. 
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cases when they are evaluated on the equipotential surface it- 
self]- Anyway, the internal charges and any other observables 
not defined on the equipotential surface, are calculated in both 
scenarios with the same set of dj coefficients. 

From the discussion above, we see that when we have a set 
of free conductors, the simplest equipotential surface that we 
can define is the one lying at infinity with zero potential, which 
is equivalent for most of the purposes to consider a cavity of a 
grounded external conductor in which all the dimensions of the 
cavity tend to infinity. 

Further, we shall see that the linear relation between charges 
and potentials in Eq. Q is gauge invariant by shifting the po- 
tential throughout the space as <p' — > <p + <po with ipo being a 
non-zero constant. This gauge transformation must keep all ob- 
servables unaltered, in particular the charge on each surface 
of the conductors. Writing Eq. ((4]) in matrix form and using Eq. 
(|10|) we have 



Q' = C E 



b ) — Ce4> = Q 



(11) 



where we used the fact that <f>o is an eigenvector of C e with null 
eigenvalue. This gauge invariance says that there is an infinite 
number of solutions (sets of potentials) for the linear equations 
(U) with given values of the charges, this fact is related in turn 
with the non-invertibility of C e . In other words, gauge invariance 
is related with the existence of an eigenvector with null eigenvalue 
which is also equivalent to the non-invertibility. On the other 
hand, the singularity of a matrix is also related with the linear 
dependence of the column (or row) vectors that constitute the 
matrix, this lack of independence in the case of C e is manifested 
in the fact that no all charges can be varied independently as 
can be seen from the expression 



Qint — 



1N+1 



(12) 



where Qint is the total charge of the internal conductors while 
Qn+i is the charge accumulated on the surface of the cavity of 
the external conductoiQ. Further, the linear dependence of the 
e-matrix can be visualized by observing that it has the same 
degrees of freedom as the r-matrix. This fact induces us to find 
expressions involving the r-matrix only. For this, we can rewrite 
Eq. @ by following the procedure that leads to Eq. 11391) 



Qk = Ckm (<Pn 



(13) 



these relations are valid for k = 1, . . . , N + 1. However, since Eq. 
(|12|) shows that Qjv+i is not independent, we can restrict them 
to k = 1, . . . , N. Rewriting Eq. (|13[) in matrix form with this 
restriction we get 



Q = CV; V=(y u Va,...,Vir), Vi = Vi -ip N+1 



(14) 



this relation is written in terms of voltages instead of potentials, 
so it is clearly gauge invariant. Further, the relation is invertible 
because the r-matrix C is positive-definite. It worths emphasiz- 
ing that all expressions obtained from now on in terms of voltages 
and the r-matrix, are valid only if the voltages are taken with 
respect to the vsjv+i potential. 



f Of course the potential on the equipotential surface is the same in 
both cases by definition. 

*This can be shown from Gauss's law or directly from the formal- 
ism presented here (see Ref. 1271 ). If the equipotential surface is a 
geometrical place in the vacuum, Eq. J 12H must be interpreted as a 
numerical equality between the total internal charge and the quantity 
on the right-hand side of Eq. J4} with i = jV + 1. 



5 Physical implications of the posi- 
tivity of the matrix 

To facilitate the derivation and interpretation of the results let 
us define the following quantities 



ko 



where ko is a constant defined such that cy are dimensionless. 
From these definitions Eq. Q could be rewritten in the form 



JV+l 



(15) 



the dimensionless Cij factors contain the same information as 
dj. Similarly, <E>i are quantities with dimension of potential 
but with the physical information of the charges Qi (it is like 
a "natural unit" for the charge). The aim of settle the charges 
and potentials with the same dimension is to interpret Eq. (|15|l 
as a linear transformation in the configuration space $ iv+1 in 
which each axis has dimensions of potential. This space would 
be isomorphic to R N+1 if we define an inner product of the form 

JV+l 



where we have taken into account that this is a real vector space. 
The capacitance matrix is hermitian (real and symmetric) with 
respect to this inner product. Now let us take two sets of charges 
and potentials and {<!>', <^'} over the same configuration 

of conductors. Doing the inner product (<&', <f>), using Eq. (|15[1 
and taking into account the hermiticity of c e , we have 



,c e </>) = (0',*) = (*,0') 



so that 



which is the reciprocity theorem. From this point of view, this 
theorem is a manifestation of the hermiticity of the e-matrix. 
Of course, we can define a potential space §> N , in which the 
N internal charges and N voltages form N— dimensional vector 
arrangements and the r-matrix acts as an hermitian operator. In 
this space the reciprocity theorem acquires the form 

(*',V) = (*,V) 

where in this case 4>' and 4> refer to configurations of the inter- 
nal charges only. Now we shall rewrite the electrostatic internal 
energy U of the system given by Eq. (JSJ in our new language 



1 

U=-, 



> 



U/k 



(16) 



the inequality comes from the positivity of the e-matrix. This 
expression is gauge invariant and can be written in terms of the 
r-matrix and voltages (see appendix [BjFl as follows 

u=i(V,cV) > 



» There is a subtlety with the concept of internal energy. The value 
of an energy is not gauge invariant, but the internal energy is indeed 
a difference of energies between an initial and a final configuration (or 
a work to ensemble a given system) this value should then be gauge 
invariant. 
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Because c is positive-definite, a zero energy is obtained only with 
V = 0. The only configurations with zero energy are the ones 
with all potentials equaO- Hence, for any geometry of the set of 
conductors and for any configuration of charges and potentials on 
them, the external agent that ensembles it, makes a net work on 
the system. There is no configuration in which the system makes 
a net work on the external agent. Note that all the analysis 
above is consistent with the features coming from the equivalent 
equation 

u= \ I E 2 dV 

where E is the electric field generated by the configuration 
throughout the volume Vs T - 

Let us interpret the eigenvalue equation of c. It reads 



cV 



(k) 



A fe V 



(k) 



(k) 



(k) 



we use superscripts to label a given eigenvector and subscripts 
to label a given component of a fixed eigenvector, if there is 
a set {i} of n indices such that all the A^'s are identical, this 
eigenvalue is n— fold degenerate. Each eigenvector V^' means a 
configuration of voltages for which each internal charge is 
related with its corresponding voltage by the same constant 
of proportionality A^. Now, since the eigenvalues are positive, 

(k) (k) 

each internal charge &\ and its corresponding voltage !/■ have 
the same sign. 

Let us construct a complete orthonormal set of real dimen- 
sionless eigenvectors uS k ' of c associated with the eigenvalues 
Xk- We show in appendix [Bl Eq. ()43p that the internal energy 
associated with a set of voltages described by the vector V can 
be written in terms of those eigenvectors and eigenvalues 



1 N 

i£A„|(u^,V 



The set |u'"'| defines principal axes in the potential space 
<£> N , and ^u*-™', is the projection of the vector V along with 

the principal axis u'™' . If the configuration of voltages in the sys- 
tem is of the form V' fc ' = Vou' fc ' (i.e. if the vector V is parallel 
to a principal axis) we fined] 



M = ^kV,f = ^A fc ||V 



(fe) 



so the eigenvalue is proportional to the internal energy associated 
with a set of voltages that forms the corresponding normalized 
eigenvector of the r-matrix. 

6 Minimization of the internal en- 
ergy 

Let us find the configuration V of voltages that minimizes the 
internal energy with the constraint that the total internal charge 



^This is in turn related with the fact that the null eigenvalue of the 
e-matrix is non-degenerate. If a degeneration of the null eigenvalue 
were present, we would have at least one eigenvector associated with 
the zero eigenvalue and linearly independent of <f>Q. The existence 
of this eigenvector would imply the existence of a configuration of 
different potentials with a null value of the internal energy. 

II Since uw are dimensionless, Vq has units of potential. Note that 
when V is parallel to a principal axis (i.e. becomes an eigenvector of 
c) , all obser vables become simpler as in the case of the axis of rotation 
in the rigid body motion. 



Qint is a constant Qo- Since Qint = —Qn+i and taking into 
account that Eq. (|13[) is also valid for k = N + 1, we have 



.V 



Qint — — y CN+l,jVj — Qo 
3=1 

the function Z (V) that defines the constraint is 



Z (V) = - £ C N+llj Vj - Qo = 

3=1 

from the Lagrange's multipliers method we have 

dU dZ „ . , 
oVi oVi 



(17) 



(18) 



(19) 



where j3 is the multiplier. Writing the internal energy as 

1 1 N N 

U=-(V,CV) = -^2^V h C kj V j (20) 

fe=U=i 

replacing Eqs. (|18U20|I into Eq. (|19(l and using the symmetry of 
the matrix, we find 



^CijVj = pC it N+i ; t = l,...,JV 

3=1 

and applying a sum over i on Eq. (|2ip 

N N N 

j=l i—l z — 1 



- £ VjCN + l,j 
3=1 



-/3Cn+i,n+i 



(21) 



(22) 



(23) 



where we have used ([5|. Substracing Eqs. (|23l I17[) and solving 
for ft we find 



/3 ; 

Eq. (|2ip can be rewritten as 



Cat 



+i,N+i 



(24) 



CV = /?v c ; = (C 1 ,n + i,C2,n+i,-,C n , n+1 ) (25) 



For a given /?, the solution of Eq. (|25[) is unique because the 
r-matrix C is invertible. It is easy to check that V = (V a , V a ), 
is a solution of Eq. (I25|l . inserting this solution in Eq. (125p . we 
get 

N 

VaY,°i3 = /9C»,iv+i ,i=l,...,N 

3 = 1 

-V a d, N+ i = pCi, N+ i , i = 1, . . . , N (26) 
where we have used ((SJ). From (|26l) we have 

Va = -P (27) 

Thus, the configuration of iV voltages that minimizes the energy 
with a fixed value of Qint, is given by 



V 1 =(-/?,-/?, -..,-#) 



(28) 



This kind of solution for V is expected because the configura- 
tion of minimal energy is obtained when we connect all the in- 
ternal conductors among them by conducting wires, this proce- 
dure clearly keeps Qint constant and equates internal potentials. 
Since all potentials of the interior conductors are the same, we 
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can define a single voltage between the external conductor and 
the internal ones, this voltage is \(3\. Since Qint = — Qjv+i = Qo 
we can figure out the system as equivalent to a system consisting 
of two conductors with charges ±Qo and voltage \/3\. Thus, we 
are led naturally to an equivalent capacitance for this system of 
potentials and charges 



|Qo| = C eq |/3| C eq = 



= Cna 



1,JV+1 



where we have used Eq. (I24[) . It can be checked that the internal 
energy U for the configuration described by (|28l) is 



1 2 1 2 

U = -CjV+l,JV + lV a = -Ceqfi 



(29) 



as expected. A brief comment with respect to the e-matrix is 
in order. This matrix has no additional degrees of freedom with 
respect to the r-matrix so that we can formally write all results 
in terms of the elements dj of the r-matrix. Notwithstanding, 
the extended elements could be useful for explicit calculations. 
Assume for instance that for the problem in the present section, 
we want to calculate the total internal charge for a given voltage 
of the system, and the equivalent capacitance. These calculations 
can be done with the following expressions 



N 



N N 



Qint — ^ ' Qk — Va 'S S ' Cij — C eq V a 
i — 1 i=l j = l 

N N N 



c e 



Cjv+i.jv+i = — Ci t M+i — Cij 



t=l 2=1 



Therefore, in terms of the elements of the r-matrix all the 
N(N + l)/2 coefficients must be evaluated with Eq. ([3]), to 
calculate Qint and C eq . In contrast, by using the e-matrix, only 
the Cn+i,n+i coefficient should be calculated through Eq. ([3]) 
to find such observables. This difference becomes more impor- 
tant as N increases. Similar advantages appear in more general 
contexts (see appendix A of Ref. [27]). 

7 The case of a chain of embedded 
conductors 

Section [3] shows that the e-matrix is singular positive and the 
r-matrix is positive-definite. The first fact was independent of 
the connectivity of Vs T . In contrast, the second fact was de- 
rived from the statement that all non-diagonal dj elements were 
strictly negative. However, it is shown in appendix \C\ that in the 
case of a chain of embedded conductors (see Fig. [2]), some of 
the non-diagonal elements dj are null because the volume Vs T 
is disconnected. Then we should check whether the r-matrix is 
still positive-definite for the chain of embedded conductors. 

Appendix [Cl shows that elements of the form d t i±i are non- 
zero in general. Appealing to an argument analogous to the one 
presented in Sec. Owe can show that for a well-behaved geometry 
of our embedded conductors, Cu > and d,i±l < 0, while the 
remaining elements vanislQ. With these properties and the fact 
that the sum of elements of each row of the e-matrix is null, we 
see that the sum of elements in each row of the r-matrix is null 
except for the N— row, for which the sum is positive. Therefore, 
the r-matrix of a chain of embedded conductors satisfies the con- 
ditions of theorem C in appendix [X] Consequently, for a chain 



of embedded conductors, the r-matrix continues being positive- 
definite, and the e-matrix is still singular positive. Combining 
these facts with theorem B of appendix [Aj we obtain that the 
null eigenvalue of the e-matrix is non-degenerat^ll It is again 
consistent with the fact that the only configuration of null inter- 
nal energy is the one associated with all conductors at the same 
potential. 



8 Conclusions 

We have studied an electrostatic system consisting of a set of 
iV conductors with an equipotential surface that encloses them. 
The associated matrix of capacitance has dimensions (N + 1) x 
(N + 1) (extended or e-matrix) even if the equipotential surface 
goes to infinity. It is usual in the literature to work with the 
matrix of dimension N x N (restricted or r-matrix), this practice 
is correct only if the voltage of the conductors is taken with 
respect to the potential of the equipotential surface. We prove 
that the e-matrix is a real positive and singular matrix, this 
is consistent with the fact that gauge invariance requires the 
existence of a null non-degenerate eigenvalue of this matrix. The 
r-matrix is a real positive-definite matrix so all its eigenvalues are 
positive. 

By constructing a "potential space" with inner product, we 
can derive some results such as the reciprocity theorem and the 
non negativity of the electrostatic internal energy of the system 
from another point of view. The eigenvectors of the r-matrix 
correspond to the sets of voltages for which such voltages and 
their associated charges on each internal conductor are related 
with the same constant of proportionality, the positivity of the 
eigenvalues guarantees that those charges and voltages are of 
the same sign. In addition, a given eigenvalue is proportional to 
the internal energy associated with the set of voltages generated 
by its corresponding eigenvector. Moreover, a complete set of 
orthonormal eigenvectors of the r-matrix defines principal axes 
in the "potential space" . 

The problem of the minimization of the internal energy is stud- 
ied under the constraint of constant value of the total internal 
charge. In this case we can define an equivalent capacitance for 
any number of internal conductors. From this problem, we real- 
ized that although the e-matrix has the same degrees of freedom 
as the r-matrix, such extension could lead to great simplifications 
of some practical calculations. 

Further, systems of successive embbeded conductors are ana- 
lyzed showing that some coefficients of capacitance are null for 
this system, allowing an important simplification for practical 
calculations. This fact is related with connexity properties of 
the volume in which Laplace's equation is considered. Moreover, 
we prove that for these configurations of embedded conductors 
the e-matrix is still positive singular with a non-degenerate eigen- 
value and the r-matrix is positive-definite. 

Finally, the properties of the matrix of capacitance shown 
here, can be useful for either a formal understanding or prac- 
tical calculations in electromagnetism. It worths observing the 
similarity in structure between the matrix of capacitance and the 
inertia tensor. 



*Of course if i = 1 (or i = N + I) the element Ci^—i (or C^i+i) 
does not exist. For a given i, at least one of them exists. 



TNote that in this case, theorem A of appendix [A] cannot be used 
to establish the non-degeneration of the null eigenvalue. 
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A Some special types of matrices 

This appendix concerns the study of a special type of matrices. 
Let define E^ as the sum of the elements on the k— row of a given 
matrix. We shall make the following 

Definition: A sp-matrix, is a square real matrix of finite 
dimension in which c m k = Ck m < for k 7^ m, and in which 
Efe is non-negative for all k. We denote with a single prime {k 1 } 
the set of indices for which E fc / > and with double prime {k"} 
the set of indices for which E fe » = 0. If no prime is used, either 
situation could happen. 

Theorem A: If C is a sp-matrix, then C is a positive matrix 
with respect to the usual complex inner product. O If {k 1 } is 
empty, the matrix is singular and vectors of the form = 
(n, . . . , n) are eigenvectors of C with null eigenvalue. Further, 
the null eigenvalue is non-degenerate if all elements of the matrix 
are non-null. © If {A:"} is empty, the matrix is positive-definite. 

Proof: We should prove that 



V T CV > 



(30) 



for an arbitrary vector V, and we should look under what condi- 
tions exists at least one non-zero vector V for which this bilinear 
expression is null. Rewriting V = N + £M with N, M being real 
vector arrangements, and using the symmetry of C, the bilinear 
form in Eq. (O becomes N T CN + M T CM. Therefore, it suf- 
fices to prove the positivity (or non negativity) of the bilinear 
form with real vector arrangements. Let N be a non-zero real 
vector, the associated bilinear form is 



N T CN = 



UkCkkiik + Ckm nkrim 

k k m^k 

for the remaining of this appendix, we assume that indices 
labeled with different symbols are strictly different. We 

rewrite the bilinear form as 



2 



^2{c kk r 

k 

+ 2 E c km[nl + n 2 m - (rife - n m ) 2 ]} 



N T CN 



E 



'kk + g E Ckm n * + \ E Ckmn rn \ 
m J m J 



- ^ ^ Ckm (rik — n m ) 2 



Now, since E fc = c kk + J2 m Ckm - °> we nave 



1^ ^ 1 



( k n 



(31) 



(32) 



2^ - 2- 

m in 

it is convenient to separate the sets {k'} and {k"} in Eq. (132[) 

c k"k" + 2 5^ Cfc " m = ~2 5^ Cfc " m (^) 
in m 

c fc'fc' + ^E Cfe '™ > _ ^E Cfc ' m ( 34 ) 



we examine first the case in which {k'} is empty. In that case 
there are no equations of the type (|34[1 , and all indices accomplish 
the equation (|33p . Using (|33|l in Eq. (|3ip and the fact that 
Cfem < 0, we find 



= E 



I -2~E Cfcm Uk + 2 E Ckmn ™ \ 
v L m J m ) 



4EE 



\Ckm \ (rik — n m ) 



using the symmetry of the matrix and taking into account that 
k, m are dumb indices, the first two terms on the right-hand side 
vanish and we find 



(35) 



Equation (|35|) shows that the bilinear form is always non- 
negative and that N T CN = for non-zero vector arrangements 
of the form = (n,n, . . . ,n). Consequently, the matrix is sin- 
gular positive. We can check that No is an eigenvector of C, with 
null eigenvalue. If all elements Ck m are non-null, Eq. (135 p shows 
that this is the only linearly independent solution, so that the 
zero eigenvalue is non-degenerate. 

Now we examine the case in which {k"} is empty, so there 
are no equations of the type (|33p , and all indices accomplish the 
equation (|34|l . Replacing (|34p into Eq. (|31|) . using the symmetry 
of C, the fact that Ck m < 0, and that N/Owe find 



N T CN > i Y, E l c *™l ( Uk ~ n '") 2 ^ 



(36) 



and the bilinear form becomes positive if and only if N 7^ 0. 
Hence, the sp-matrix is positive definite when {k"} is empty. 
QED. 

Theorem B: Let C be a matrix of dimension (jV + 1) x 
(TV + 1), such that E; = for all rows. This matrix has eigen- 
vectors of the form Njf = (n, . . . , n) associated with a null eigen- 
value. Let C r be the N x N submatrix of C consisting of the 
elements dj of C with i, j — 1, . . . , N . If C r has no null eigen- 
value^, the null eigenvalue of C is non-degenerate. 

Proof: The condition Ei = for i = 1, . . . , N + 1 gives 



JV+l 

E c - k = 



i = L 



,N + 1 

Eigenvectors of C with null eigenvalues must give 

; i= 1,...,N + 1 



JV+l 



(37) 



(38) 



assuming Uk — n for all k and using condition (I37|l . Eq. 
is satisfied. Thus, Njf is an eigenvector associated with a null 
eigenvalue. From the condition (|37p we also find 



JV+l 

CikUk 

fc=i 



= CikTik + Ci,JV+lJlJV+l 



JV+l 

C ik n k 



k=l 
N 

CikUk - 

JV 

^CikK-nw+i); i = l,..,N + l (39) 



- 2^ Cik 1 " \ i 1 



If C r is a normal matrix (or if it can be brought to the canonical 
form), it is equivalent to say that C r is non-singular. 
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replacing ([39]) in 



the latter becomes 



njv+i) = 0; i = l,...,N + l 



(40) 



in particular Eq. (I40[) holds for £ 
Eq. (gOJ) becomes 



, N. With this restriction 



C r V = ; V = (ni - njv+i, n 2 - njv+i, njv - njv+i) (41) 

since C r has no null eigenvalues, the only solution for Eq. (|4ip 
is V = 0. Hence the only type of solutions for N are of the form 
N = (njv+i, njv+i, tijv+i) which are all linearly dependent. 
Hence, the null eigenvalue is non-degenerate. It is immediate 
that these solutions satisfy Eq. (|40[) for i = N + 1 as well. Note 
that C is not neccesarily symmetric or real. QED. 

Theorem C: Let C be a N x N sp-matrix such that {k"} = 
{1, 2, . . . , N — 1} and the terms 



s+1 — Ci+1 



2 , <-Z,l 



-1 = Ci-1A 



i = 2. 



are non-zero, while the remaining non-diagonal terms vanish. 
Then C is positive-definite. 

Proof: Assume C as singular and arrive to a contradiction. 
Replacing Eqs. ([33] [34]) in Eq. J3TJ we obtain 



{—— 2 i i. 2 I 

m J m J 



i££ Cfcr 



2 

n m ) 



such a replacement also shows that the equality holds if njv = 0, 
while the strict inequality holds if njv 7^ 0. Using the symmetry 
of the matrix and the fact that ctm < we have 



1 



Since C is singular, a non-trivial solution must exist for the bi- 
linear form to be null. For this, the equality must hold in this 
relation, therefore njv = 0. The term on the right written in 
terms of the non-zero elements of the matrix yields 



N T CN 



C12 

JV-1 

+ £ 



n 2 f + £ 



|Cfc,fc+l| , ,2 
(71* - Tlfc+i) 



|Cfc,fc-i| 



rik-i 



) 2 + 



-n N _ 



for this expression to be zero each term in these sums must be 
zero. Since all matrix elements involved in this expression are 
non-zero, the last sum says that njv-i = 0, while the other sums 
say that 711=712 = ... = tijv-i- Since njv was already zero, it 
shows that the only solution is the trivial one, contradicting the 
singularity of the matrix. QED. 

B Some properties of the internal 
energy 

The equation (|16p for the internal energy u in "natural units" 
can be written in terms of voltages instead of potentials with the 
r-matrix. By using the fact that (j>o is an eigenvector of c e with 
null eigenvalue, and the hemiticity of c e , Eq. (|16|) becomes 



2 it 
2 it 



(c e 



-&>)) = ( c <=V»i (<t>-<t>o)) 

b o) . (<t> ~ 4>o)) = {<t> ~ 00: C = (<t> ~ <t>o)) 



defining a vector arrangement of N + 1 voltages V e we find 

2u = (V e ,c e V e ); Vj = {V 1 ,...,V N ,V N+1 ), Vi = <pi- <po 

writing this bilinear form explicitly and expanding the sums we 
find 

. JV+l JV+l 

u = 2 £ £ CijViVj 
i=i 3=1 

^ JV JV 

= 2 ^ ^ + l/iY + lifjv + 1 

i=l 3=1 

* 1 
Kn+i = Ci,jv+iVi + -cjv+i,jv+i Vat+1 

i=i 

choosing tpo = <£>jv+i we find Vjv+i = 0, hence 

JV JV 

u = 2^^ Cl]VlV ^ 2 (V ' CV) 
i=i j=i 

V T = (Vi,...,Vat) , Vi = y>i — v^iv+i 

thus the internal energy can be written in terms of the r-matrix 
and the voltages in a simple way as long as the latter are defined 
with respect to the external potential ipN+i- Since V is gauge 
invariant u also is. 

On the other hand, remembering that we can always construct 
a complete orthonormal set of real dimensionless eigenvectors 
Ue of the e- matrix associated with the eigenvalues X%, we can 
write the internal energy associated with a configuration <f> of 
potentials in terms of these eigenvalues and eigenvectors. Since 
the eigenvectors form a basis we can express cj> as a linear com- 
bination of them 



^ b m ui m > ; b m 



and Eq. (|16l) becomes 



JV+l 



JV+l 



2u = {4>,c e (f>) = ^2 6 m u^ m) ,c c ^2 fe„ui n) 

\m=l 

JV+l JV+l 
m=l n=l 

JV+l JV+l JV+l 

= y ] y ] XnbmbnSmn — ^ ] X n b n 
m—1 n — 1 n—1 

1 N + 1 



(42) 



Since the eigenvectors u' fc ' and the matrix c e are dimensionless, 
the eigenvalues AJ: also are. It is straightforward to write this 
expression in terms of the r-matrix and the voltages with respect 

to IfN+l 

1 N 2 



where X n ,u^ are eigenvalues and eigenvectors of the r-matrix 

C Some properties of chains of em- 
bedded conductors 

Let us study a set of N + 1 conductors which are successivily 
embedded. We label them fc = 1, . . . , N + 1 from the inner to 
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Figure 2: A chain of embedded conductors with N = 3. 
The surfaces S2,Sa, have an inner and an outer part. 

the outer. Observe that the surface Sk for each conductor with 
k — 2, . . . ,N has an inner and an outer part, but for Sn+i we 
only define an inner part and for Si we only define an outer part 
(see Fig. [2}. In addition, we define Vk with k = 1, . . . , N as the 
volume formed by the points exterior to the conductor k and 
interior to the cavity of the conductor k + 1 that contains the 
conductor k. Let us examine the non-diagonal elements Ckm 
assuming from now on that k < m. 

From Eq. ([T]) we see that if m — k = 1 then f m (Sk) = and 
fm (Sk+i) = 1 because Sk+i ~ S m , the volume Vk is precisely 
delimited by the outer part of the surface Sk and the inner part 
of the surface Sk+i; thus / m has a non-trivial solution in Vk- 
Therefore, we have in general that V/ m 7^ in Vk and in the 
surfaces that delimite it. Thus the integral 

C km = so j V/ m ■ n fe dS (44) 

S k 

has a contribution from the outer part of Sk- Now, if Vk-i 
exists (i.e. if k > 1), and taking into account that f m (Sk-i) 
/,, (Sk) = 0, the uniqueness theorem says that the only solution 
in Vk-i is fm = and hence V/ m = in this volume and in the 
surfaces that delimite such a volum^E Thus the integral surface 
in (|44|l has no contributions from the inner part of Sk- 

Now, if m — k > 2 we see that f m (Sk) = fm (Sk+i) = 0, then 
the only solution in Vk is fm = V/ m = in this volume and in the 
surfaces that delimite such a volume. Thus the integral surface in 
(|44[) has no contributions from the inner part of Sk ■ On the other 
hand, if Vk-i exists (k > 1), and since f m (Sk-i) = fm (Sk) = 
we see once again that f m = V/ m = in the volume Vk-i 
and in the surfaces that delimite it; so the integral (|44p has no 
contribution from the outer part of Sk either. 

From the previous discussion and appealing to the symmetry 
of the e-matrix, we conclude that Ckm = for \m — k\ > 2. In 
addition, when \m — k\ = 1, the surface integral (|44[) receives 
contribution only from the outer part of Sk- Notice that the 
previous behavior has to do with the fact that the total volume 
Vs T consists of several disjoint (and so disconnected) regions and 
that \k — m\ > 2 indicates that these labels are always associated 
with disjoint volumes. In the last discussion we have not included 
the possibility that the most interior conductor has a cavity. 
Since it would be an empty cavity, the surface and volume of 
this cavity do not contribute to the calculation of any coefficient 
of capacitance (see Ref. |27|). 

From the results above, we see that for N + 1 successivily 



§ Remember that the surfaces are slightly different from the surfaces 
of the conductors for the gradient to be well-defined. 



embedded conductors with N > 2, we have 

Cn — — (Ci,i_i + Ci,i+i) ; i = 2, . . . , N 
Cn = — C12 ; C"jv+i,jv+i = — Cn,n+i 

How many degrees of freedom do we have for the e-matrix?. 

D Suggested Problems 

For checking the comprehension of the present formulation and 
its advantages, we give some general suggestions for the reader. 

1. Show all the properties stated here for the r-matrix and 
e-matrix with specific examples. 

2. Look for differences and similarities between the matrix of 
capacitance in electrostatics and the inertia tensor in me- 
chanics, from the physical and mathematical point of view. 

3. From (<j>,c<f>) = k we find (u, cu) = 1 with u = 4>/Vk- 
This defines the equation of an ellipsoid, describe how to 
find the length of the axes of the ellipsoid in the & N and 
3 >JV+1 spaces for the r-matrix and the e-matrix respectively. 
Describe the principal axes in these "potential spaces" . 

4. By setting du/difi — 0, prove that in the absence of con- 
straints, the only local minimum of the internal energy is 
given by sets of the type (f>o- 

5. Prove Eq. I|29p for the minimal internal energy under the 
constraint of constant internal charge. 

6. Let {a, 6} be two positive numbers. Consider the 4x4 
matrix given by 

„ ( aB 2x2 2x2 \ „ _ / 1 -1 \ 

C4X4= (, 2x2 6B 2X2 ) ; B2X2 = (, -1 1 ) 

this is a sp-matrix in which the sum of elements in each row 
is zero. Further, A — is a two-fold degenerate eigenvalue 
of C. Can C be a matrix of capacitance associated with a 
given electrostatic set of conductors?. 

7. Look up for more applications of singular positive and 
positive-definite matrices in different contexts of Physics. 
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